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Abstract 


The nucleotide sequences of the internal transcribed spacers (ITS 1 and ITS 2) and 5.8S coding region of nuclear 
ribosomal DNA, as well as the non-coding trnL-trnF spacer regions of the chloroplast DNA, were determined and 
analyzed to estimate the phylogenetic position of Takhtajania perrieri (Capuron) Baranova & J.-F. Leroy within the 
Winteraceae. Using representatives of each genus of Canellaceae as outgroups ( Canella , Capsicodendron , Cinnamoden- 
dron , Cinnamosma , Pleodendron , and Warburgia), both maximum parsimony and maximum likelihood analyses place 
Takhtajania in a basal position sister to the remainder of the Winteraceae. Although the overall topology within the 
Winteraceae was mostly congruent between nuclear and chloroplast data sets, the trnL-trnV data resulted in lower 
support values in comparison to the ITS data, and failed to resolve basal relationships in the family, yielding alternative 
equally parsimonious solutions. The combined nuelear/chloroplast data set resulted in a single tree identical to that 
generated by nuclear data alone, but demonstrated strengthened support for a basal branch leading to Takhtajania , as 
well as for the position of Tasmannia as the next branch within the family. Potentially conflicting signals from nuclear 
and chloroplast data indicate that further taxon sampling or additional sequence data may be required to infer infra- 
familial phylogenetic relationships for Canellaceae. 

Key words: Canellaceae, combined data, cpDNA, ITS, Likelihood, phylogeny, Takhtajania , W interaceae. 


Known until recently from only a single collec¬ 
tion made in 1909, Takhtajania perrieri has long 
been an enigma among the basal angiosperm family 

Winteraceae (Leroy, 1978). Its rediscovery in 1994 
(Schatz et al., 1998) has provided the opportunity 
to reevaluate those morphological features anoma¬ 
lous within the family, as well as to explore the 
phylogenetic relationship of this isolated, relictual 
Malagasy endemic. Based upon morphological 
characters, two alternative hypotheses have been 
proposed regarding the position of Takhtajania 
within the Winteraceae: (1) Takhtajania is a mem¬ 
ber of a clade also comprising Pseudowinter a, Bel - 
liolum , Bubhia , Exospermum , and Zygogynum (the 


latter 4 genera sometimes combined into a single 
genus Zygogynum s.l.) (Vink, 1988; Endress et al., 
2000 this issue); or (2) Takhtajania represents the 
sister taxon to the remainder of the Winteraceae 
(Leroy, 1978; Vink, 1988). The aim of the present 
study is to infer the phylogenetic position of Takh¬ 
tajania based upon analyses of molecular sequence 
data. 

A previous molecular phylogenetic study (Suh et 
al., 1993) utilized the internal transcribed spacer 
(ITS) region of the nuclear ribosomal DNA (nrDNA) 
to address questions of generic relationships in the 
Winteraceae. Tissue of Takhtajania , however, was 
unavailable. In addition to presenting a phylogeny 
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for the family solely through ingroup analysis, Suh 
et al. (1993) identified a single gene duplication 
event in the Zygogynum s.l. clade, a possible mo¬ 
lecular marker for the inclusion (or exclusion) of 
Takhtajania within the Belliolum/Bubbia/Exosper- 
mum/Zygogynum assemblage. 

If an ITS duplication event either was not iden¬ 
tified in Takhtajania or occurred independently, 
then outgroup analysis would become critical for 
investigating its relationship to the remainder of the 
family. Both morphological and molecular (r6cL) 
cladistic studies, which examined relationships 
among basal angiosperms, have identified Canel- 
laceae (represented in these studies by Canella 
alone) as a putative sister taxon to Winteraceae 
(Donoghue & Doyle, 1989; Qiu et al., 1993; Chase 
et al., 1993). Canella also appeared to be closely 
related to Winteraceae on the basis of secondary 
metabolites (Gottlieb et al., 1989). Canellaceae are 
comprised of six genera found in tropical Africa 
(Warburgia ), Madagascar ( Cinnamosma ), and trop¬ 
ical America ( Capsicodendron and Cinnamoden- 
dron in South America, Canella and Pleodendron 
in the Caribbean) (Cronquist, 1981). Although Suh 
et al. (1993) had encountered difficulties in the 
alignment of ITS sequences between Canella and 
Winteraceae, it was decided to reattempt alignment 
after sequencing additional representatives of Ca¬ 
nellaceae. 

This paper reports on the implications of new 
ITS sequence data for Takhtajania and the five oth¬ 
er genera of Canellaceae; Canella ITS sequence 
had been previously determined by Suh et al. 
(1992). In addition, sequence data have been gen¬ 
erated for two tandem chloroplast encoded (cp- 
DNA) spacer regions, and analyzed both separately 
and in combination with the ITS sequences. The 
first spacer is a group I intron and is located be¬ 
tween the conserved trnL (UAA) 3’ exon and trnh 
(UAA) 5' exon. The second, an intergenic spacer, 
is located between the trnh (UAA) 5' exon and trnF 
(GAA) (Taberlet et al., 1991). Restriction fragment 
length polymorphisms of this molecule have been 
used to address species-level relationships within 
Mexican pines (Perez de la Rosa et al., 1995). Se¬ 
quence data generated from the trn\-Y region also 
have been utilized to resolve phylogenetic relation¬ 
ships at both the species level (Mes & Hart, 1994) 
and the generic level (Compton et al., 1998; 
McDade & Moody, 1999). Several studies (Gielly 
et al., 1996; Molvray et al., 1999) have demonstrat¬ 
ed the utility of combining the ITS and trnh-F spac¬ 
er sequences for phylogenetic reconstruction. Thus, 
a complete survey of both ITS and trnh-F regions 
for all genera of both Winteraceae and Canellaceae 


has been conducted to estimate the phylogenetic 
position of Takhtajania. 

Materials and Methods 

SPECIES SAMPLED AND SOURCES OF PLANT 
MATERIAL 

The species included in this study are presented 
in Table 1, along with voucher information, litera¬ 
ture citation, and GenBank accession numbers for 
both the ITS and trnh-F spacer sequences. This list 
includes at least one representative species from all 
genera in both the Winteraceae and Canellaceae. 
In the instances where ITS has been determined 
previously (all Winteraceae [Suh et al., 1993] ex¬ 
cept Takhtajania, plus Canella [Suh et al., 1992]), 
identical DNA extracts were used to amplify and 
sequence the trnL-F spacers. 

DNA EXTRACTION, AMPLIFICATION, AND 
SEQUENCING 

Standard CTAB methods of DNA isolation (Doyle 
& Doyle, 1987) were used to obtain total genomic 
DNA for polymerase chain reaction (PCR) ampli¬ 
fication. To generate ITS amplicons, PCR was per¬ 
formed using plant specific primer ITS5 (Suh et al., 
1993) and universal primer ITS4 (White et al., 
1990) under the conditions described in Suh et al. 
(1993). To generate the trnh-F spacer amplicons, 
PCR was performed using universal primers “c” 
and 4 T’ (Taberlet et al., 1991), also under the con¬ 
ditions described in Suh et al. (1993). 

Sequence data were generated from PEG (poly¬ 
ethylene glycol) purified double-stranded PCR 
products (Morgan & Soltis, 1993) with the dye-ter¬ 
minator cycle-sequencing protocol for the ABI 
373A Sequencer (Applied Biosystems, Inc.). The 
entire forward and reverse strand of ITS was deter¬ 
mined for each species using the amplification 
primers (see above) and primers ITS3 (White et al., 
1990) and C5.8S (Suh et al., 1993). To determine 

the sequence for both strands of the trnL-F spacers, 
the amplification primers were used as well as 
primers “d ,? and “e” (Taberlet et al., 1991). The 
resulting chromatograms were edited with Se- 
quencher version 3.1 (Gene Codes, Inc.); regions 
corresponding to the amplification primers were de¬ 
leted. The final edited consensus sequences were 
exported for alignment. 

SEQUENCE ALIGNMENT 

All sequences were aligned manually with the 
aid of Se-Al version l.Oal (Rambaut, 1996) mul¬ 
tiple sequence editor. 
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lahle 1. List of specimens used in this study with locality and voucher information, ITS literature citation for 
previously published sequences, and GenBank accession numbers for both ITS and trnh-F. 


Species 

Winteraceae 

Belliolum pancheri (Baill.) 

Vink (= Zygogynum panch¬ 
eri (Baill.) Vink) 

Buhbia comptonii (Baker f.) 
Dandy (= Zygogynum 
comptonii (Baker f.) Vink) 
Drimys winteri J.R. Forst. & G. 
Forst. 

Exosperm urn stipitatum (Baill.) 
Tiegh. ex Morot (= Zygo¬ 
gynum stipitatum Baill.) 
Pseudotvin tera axilla ris (J. B. 

Forst. & G. Forst.) Dandy 
Pseudowintera colorata (Raoul) 
Dandy 

j 

Takhtajania perrieri (Capuron) 
Baranova & J.-F. Leroy 

Tasmannia insipida R.Br. ex 

DC. 

Tasmannia lanceolata (Poir.) 

A.C. Smith 

Zygogynum acsmithii Vink 

Zygogynum balansae Tiegh. 

(= Zygogynum pom ifer urn 
Baill. subsp. balansae 
(Tiegh.) Vink) 

Zygogynum bicolor Tiegh. 

Canellaceae 

Canella winterana (L.) Gaertn. 
Capsicodendron din is ii 
(Schwacke) Occhioni 
Cinnamodendron ekmanii Sleu- 
mer 

Cinnarnosrna madagascariensis 
Danguy 

Pleodendron m acro nth u rn 
(Baill.) Tiegh. 

Warburgia salutaris (G. Bertol.) 
Chiov. 


Voucher 


New Caledonia, Ml. koghis, 

LBT 205, NO 

New Caledonia, Prokomeo, NW 

of Canala, LBT 204 , NO 

South America, Chile, (Berke¬ 
ley B. G.) 45.307, NO 

New Caledonia, Mt. Pani£, LBT 

202 , NO 

New Zealand, N. Island, Aka- 
tarawa, LBT 300, NO 
New Zealand, N. Island, Aka- 
tarawa, LBT 301 , NO 

Madagascar, Anjahanaribe-Sud 
Special Reserve, Rakotomal- 
aza et ai 1342 , MO 
Australia, Queensland, Ather¬ 
ton, LBT 108, NO 
Australia, (Berkeley B. G.) 

60.0052, NO 

New Caledonia, Mts. near Lac 
Enhuit, LBT 201, NO 

New Caledonia, Ozurnae Mts., 

LBT 203, NO 


New Caledonia, Plateau de 
Dogny, LBT 200, NO 


South America, LBT 124 , NO 
South America, Butzke et ai 

11521 , US 

Dominican Republic, Samana 
Penn insula, Garcia & Yeloz 

6866, SD 

Madagascar, Lowry 4001, MO 
Puerto Rico, Axelrod 10783, 

MO 

South Africa, Mpumalanga, 

Goldblatt 11314, MO 


Genbank 


ITS citation 

Genbank ITS 

/niL-F 

Suh et al., 1993 

AY004119 

AY004135 

Suh et al., 1993 

AY004123 

AY004140 

Suh et al., 1993 

AY004126 

AY004143 

Suh et al., 1993 

AY004121 

AY004138 

Suh et ah, 1993 

AY004124 

AY004141 

Suh et al., 1993 

AY004125 

AY004142 

This paper 

AY004129 

AY004146 

Suh et al., 1993 

AY004127 

AY004144 

Suh et al., 1993 

AY004128 

AY004145 

Suh et ah, 1993 

AY004122 

AY004139 

Suh et al., 1993 

AY004120 

AY004137 

Suh et al., 1993 

AY004118 

AY004135 

Suh et al., 1992 

L03844 

AY004152 

This paper 

AY004132 

AY004149 

This paper 

AY(X)4133 

AY004150 

Hi is paper 

AY004131 

AY004148 

This paper 

AY004134 

AY004151 

This paper 

AY004130 

AY004147 


Nuclear ITS . The alignment of Winteraceae 
ITS sequences of Suh et al. (1993) was used as a 
guide to easily incorporate the Takhtajania ITS se¬ 
quence. The paralogous copies (GenBank 
AY004111—AY004117) of the ITS region as deter¬ 
mined by Buckler et al. (1997) were not included 
in the alignment. 


Alignment of ITS sequences within Canellaceae 
also was easily determined by eye. Alignment of 
Canellaceae with Winteraceae sequences into a 
single data set required additional information to 
generate a plausible alignment. One region of ITS 
1 and another in ITS 2 were unalignable across 
both families, though these regions were alignable 
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among families. These regions were treated as non¬ 
overlapping insertion/deletion (indel) events in the 
final data set. Inferred indel regions were included 
in all phylogenetic analyses with the resulting gaps 
treated as missing. 

Putative secondary structures of ITS 2 were de¬ 
termined for each taxon employing the minimum 
free-energy program mFOLD (Zuker, 1989) using 
constraints previously described by Hershkovitz 
and Zimmer (1996). In their study, the putative sec¬ 
ondary structure of Canella ITS 2 was determined, 
and this structure was used as a guide to determine 
structures for all remaining Canellaceae. 

Chloroplast tm L-F. The trn L-F spacer regions 
were easily aligned by eye within and between both 
families, with the resulting indels treated identi¬ 
cally to those in the ITS data set. 

Combined data set. The aligned ITS and trn L- 
F spacer sequences were then combined into a sin¬ 
gle data set. Two partitions were defined (ITS and 
trnL-F) and subjected to the data partition-homo¬ 
geneity test (Farris et al., 1995). A thousand rep¬ 
licates were performed and the resulting P-value 
was used to determine if using the combined data 
set for phylogenetic reconstruction would be appro¬ 
priate. 

In summary, three data sets were examined: (1) 
the ITS data set consisting of the ITS 1, 5.8S, and 
ITS 2 regions; (2) the tmL-F data set consisting of 
a group I intron, the 5'-fraL exon, a non-coding 
intergenic spacer region, and partial sequence of 
trnF; and (3) a combined data set consisting of both 

the ITS and trnL-F regions. 

PHYLOGENETIC ANALYSES 

All phylogenetic reconstruction analyses were 
conducted using version 4.0b 1 or 4.0b2 of PA UP* 
(Swofford, 1998). 

Maximum parsimony. For each data set, phy¬ 
logenetic reconstruction under maximum parsimo¬ 
ny (MP) was conducted by utilizing the Branch and 
Bound search option in PAUP* with TBR branch¬ 
swapping, MULPARS, and ACCTRAN options ac¬ 
tive. Characters were assigned equal weights at all 
nucleotide positions (Fitch, 1971). Robustness of 
cladistic linkages was evaluated with 1000 boot¬ 
strap replicates (Felsenstein, 1985; Sanderson, 
1989). Decay values (Bremer, 1988; Donoghue et 
al., 1992) for each node were calculated using the 

TOPOLOGICAL CONSTRAINTS option in PAUP*. 

Maximum likelihood. The strategy employed 
utilized the likelihood-ratio test statistic (Felsen¬ 
stein, 1981; Goldman, 1993; Yang et al., 1995) to 

determine the best model of DNA substitution for 


each of the three data sets described above. The 
models considered include the general time-revers¬ 
ible model (GTR, equals REV of Yang, 1994a), 
Hasegawa et al. (1985; denoted HKY85), Kimura 
(1980; denoted K2P), and Jukes and Cantor (1969; 
denoted JC69). Equations for calculating likeli¬ 
hoods under these DNA substitution models are 
given in Swofford et al. (1996). An iterative pro¬ 
cedure to first evaluate models and optimize model 
parameters was used on an initial set of trees gen¬ 
erated from the maximum parsimony analysis. 

For each model of DNA substitution, models that 
incorporate rate variation across sites were ex¬ 
plored. These included equal rates assumed across 
all sites, a proportion of sites assumed to be in¬ 
variable with equal rates assumed across variable 
sites (“I” [Hasegawa et al., 1985]), all sites as¬ 
sumed to follow a discrete approximation of the 
gamma distribution (“T” [Yang, 1994b]), and some 
sites assumed to be invariable with gamma distrib¬ 
uted rates at variable sites (“I + P’ [Gu et al., 

1995; Waddell & Penny, 1996]). 

For the combined data set, a set of models was 
tested which explored the possibility that the rel¬ 
ative rates between the two molecules differed sig¬ 
nificantly. This class of models assigns sites to clas¬ 
ses and then estimates the relative rate for each 
class separately (Felsenstein, 1991; Goldman & 
Yang, 1994). Sites were assigned to two classes (ITS 
and trnL-F). This model is denoted as “+ SS.” 

For each data set, a maximum likelihood heuris¬ 
tic search with TBR branch swapping was then per¬ 
formed using parameters estimated from the tree 
with the best likelihood score under the best model. 
If the resulting tree was different in topology from 
that of the original tree, the resulting tree was used 
to further optimize the model parameters. With pa¬ 
rameters re-optimized, a heuristic search of ten rep¬ 
etitions of random taxon addition and TBR branch 
swapping was performed using fully defined model 
parameters. Bootstrap methods (Felsenstein, 1985; 
Sanderson, 1989) with 1000 replicates were per¬ 
formed to estimate robustness of nodal support. 

Results 

The combined data set used in this study can be 
found at TreeBASE (http://www.herbaria.harvard, 
edu/treebase/). 

SEQUENCE ALIGNMENT 

Nuclear ITS. ITS sequence alignment proved to 
be a relatively simple task within both Winteraceae 
and Canellaceae. With the exception of two regions, 

one in ITS 1 (bp 1143—1163) and the other in ITS 
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Table 2. List of the insertion/deletion (indel) events for tlie /r/iL-F/ITS combined data set. Character name (A—X 
for trn L-F; 1-115 for ITS), position information (character number in the combined data set), indel length (bp), character 
type (I = parsimony-informative, U = uninformative, H = homoplastic when mapped on Figs. 1 and 2a), and the 
species where indels were found are noted. Bp = Belliolum pancheri , Be = Bubbia comptonii , Dw = Drimys winteri , 
Ks = Exospermum stipitatum, Pa = Pseudowintera axillaris. Pc = Pseudowintera color at a, Tp = Takhtajania perrieri, 
Ti = Tasmannia insip ida, T1 = Tasmannia lanceolata , Za = Zygogynum acsmithii, Zba = Zygogynum ba I ansae, Zbi 
= Zygogynum bicolor. Cw = Canella winter ana, Cd = Capsicodendron dinisii, Ce = Cinnamodendron ekmanii, Cm = 
Cinnamosma madagascariensis, Pm = Pleodendron macranthum , Ws = Warburgia salutaris . 


Character 

Position in combined 

data set 

Length (bp) 

1/U/H 

Clade 

ImL-F 

A 

114 

i 

I 

Winteraceae/Canellaceae 

B 

179-181 

3 

I 

W interaceae/Canellaceae 

C 

207-212 

6 

U 

Tp 

[) 

247-255 

9 

u 

Cw 

E 

252-256 

5 

1 

Winteraceae/Canellaceae 

F 

257-264 

8 

u 

Ws 

G 

282-286 

4 

u 

Cm 

H 

317 

l 

u 

Cw 

I 

319-331 

13 

u 

Cw 

J 

323-331 

9 

u 

Tp 

K (a, b, c) 

359-364 

i 

H 

a = (Ti, TI, Tp), b = (Tl, Tp), 
c = Tl 

L 

379-384 

6 

1 

Dw, Pc, Pa, Be, Bp, Es, Za, 

Zba, Zbi 

M 

385-388 

4 

U 

Cd 

N 

414—418 

5 

I 

Winteraceae/Canellaceae 

O 

586-591 

6 

H 

Hun of Adenines (see text) 

P 

590-694 

105 

1 

Wi nteraceae/Canel laceae 

Q 

633 

1 

u 

Dw 

B 

634 

1 

LI 

Tp 

s 

706 

1 

I 

Winteraceae/Canellaceae 

T 

710 

1 

I 

W i nteraceae/Canel laceae 

U 

717-718 

2 

I 

Wi nteraceae/Canel laceae 

V 

740-750 

11 

U 

Ws 

w 

772 

1 

u 

Ti 

X 

897 

1 

H 

Cd, Cm, Ws 

ITS-1 

1 

993 

1 

I 

Winteraceae/Canellaceae 

2 

997 

1 

1 

Be, Bp, Ks, Za, Zba, Zbi 

3 

998-999 

2 

1 

Cd, Ce, Cm, Pm, Ws 

4 

1000 

1 

U 

Cd 

5 

1003-1013 

11 

U 

Cd 

6 

1002 

1 

1 

Pa, Pc 

7 

1003 

1 

U 

Tp 

8 

1008 

1 

1 

W interaceae/Canellaceae 

9 

1009 

1 

u 

Cw 

10 

1027-1028 

2 

I 

W i n t e raeeae/( ^ a nel 1 aceae 

11 

1040 

1 

u 

Dw 

12 

1041 

1 

I 

Wi nteraceae/Canel laceae 

13 

1042 

1 

H 

W interaceae, Cd 

14 

1045 

1 

I 

Cd, Cm, Ws 

15 

1048-1049 

2 

H 

Canel laceae, Dw 

16 

1050-1051 

2 

1 

Winteraceae/Canellaceae 

17 

1052 

1 

1 

Cd, Cm, Ws 

18 

1053-1054 

2 

U 

Ws 

19 

1058 

1 

H 

Canel laceae, Dw 

20 

1059 

1 

H 

Cd, Dw 
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Table 2. 

Continued. 




Character 

Position in combined 

data set 

Length (bp) 

I/U/H 

Clade 

21 

1060-1061 

2 

U 

Dw 

22 

1068 

1 

u 

Cw 

23 

1076-1077 

2 

u 

Cw 

24 

1081-1082 

2 

I 

Tl, Ti 

25 

1083-1084 

2 

I 

Winteraceae/Canellaceae 

26 

1090 

1 

H 

Winteraceae, Cm, Pm 

27 

1091 

1 

H 

Ce, Cw, Pm 

28 

1094—1106 

13 

U 

Cd 

29 

1194-1097 

4 

u 

Ws 

30 

1198-1104 

7 

H 

Ce, Cw, Pm, Ws 

31 

1105-1106 

2 

U 

Pm 

32 

1112-1113 

2 

u 

Dw 

33 

1114 

1 

I 

Winteraceae, Cw 

34 

1115 

1 

I 

Winteraceae/Canellaceae 

35 

1116-1118 

3 

1 

Winteraceae, Cw 

36 

1119-1122 

4 

I 

Winteraceae/Canellaceae 

37 

1128 

1 

u 

Cw 

38 

1135 

1 

1 

Winteraceae/Canellaceae 

39 

1141-1142 

2 

H 

Winteraceae, Cd, Cw 

40 

1143-1146 

4 

I 

Winteraceae/Canellaceae 

41 

1147-1163 

17 

I 

Winteraceae/Canellaceae 

42 

1148-1149 

2 

U 

Cd 

43 

1158 

1 

I 

Cm, Ws 

44 

1176 

1 

u 

Cw 

45 

1185-1188 

4 

1 

Winteraceae/Canellaceae 

46 

1199 

1 

u 

Dw 

47 

1206-1225 

20 

1 

Winteraceae/Canellaceae 

48 

1209-1212 

4 

1 

Ce, Pm 

49 

1223-1225 

3 

I 

Ce, Pm 

50 

1234 

1 

I 

Winteraceae/Canellaceae 

51 

1244 

1 

u 

Cw 

52 

1245 

1 

I 

Winteraceae/Canellaceae 

53 

1262-1277 

16 

I 

Winteraceae/Canellaceae 

54 

1262-1263 

2 

I 

Ce, Pm 

55 

1272 

1 

1 

Cd, Cm, Ws 

56 

1273 

1 

I 

Ce, Pm 

57 

1283-1284 

2 

H 

Ce, Cw, Pm 

58 

1287 

1 

H 

Cd, Cm, Cw, Ws 

59 

1292-1293 

2 

I 

Be, Bp, Ks, Za, Zba, Zbi 

60 

1299 

1 

H 

Cm, Pm 

61 

1308-1326 

19 

I 

Winteraceae/Canellaceae 

62 

1311 

1 

I 

Bp, Es. Za, Zba, Zbi 

63 

1312-1313 

2 

H 

Be, Pa, Ti, Tl, Tp 

64 

1314-1316 

3 

U 

T P 

65 

1325 

1 

H 

Winteraceae, Ws 

66 

1326 

1 

H 

Cd, Cm, Cw 

67 

1331 

1 

U 

Cw 

5.8S 

68 

1487 

1 

U 

Cw 

ITS-2 

69 

1498 

1 

H 

Pa, Pc, Be, Es, Za, Zba, Zbi 

70 

1504 

1 

I 

Winteraceae/Canellaceae 

71 

1511-1513 

3 

1 

Winteraceae/Canellaceae 

72 

1514 

1 

H 

Cd, Ce, Cm, Pm 

73 

1523 

1 

H 

Cd. Ce, Cw, Pm, Ws 
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Table 2. 

Continued. 






Posit 

ion in combined 




Character 


data set 

Length (bp) 

l/U/H 

Clade 

74 


1524-1525 

2 

H 

Be, Bp, Cm, Ks, Ti, Tl, Za, 






Zba, Zbi 

75 


1526 

1 

I 

Be, Bp, Cm, Ks, Za, Zba, Zbi 

76 


1527 

1 

I 

Winteraceae/! Canellaceae 

77 


1534 

1 

U 

Dw 

78 


1535-1539 

5 

I 

Winteraceae/Canellaceae 

79 


1540 

1 

H 

Canellaceae, Be 

80 


1543-1544 

2 

U 

Cw 

81 


1545 

1 

I 

Winteraceae/Canellaceae 

82 


1553-1554 

2 

1 

Pa, Pc 

83 


1573-1575 

3 

I 

Winteraceae/Canellaceae 

84 


1595-1596 

2 

I 

Bp, Zb 

85 


1603 

1 

H 

Be, Bp, Ce, Cm, Ks, Pm, Ti, Tl, 






Za. Zba, Zbi 

86 


1604 

1 

1 

Be, Bp, Cm, Ks, Ti, Tl, Za, 






Zba, Zbi 

87 


1612-1618 

7 

U 

Cw 

88 


1651 

1 

U 

Tp 

89 


1652-1653 

2 

I 

W interaceae/Canellaceae 

90 


1659 

1 

I 

Winteraceae/Canellaceae 

91 


1662 

1 

U 

Pm 

92 


1663 

1 

I 

W i nteraceae/Canel laceae 

93 


1672 

1 

H 

Winteraceae, Cm 

94 


1676 

1 

I 

Cd, Cm, Ws 

95 


1683-1688 

6 

U 

Ce 

96 


1690-1691 

1 

1 

Pa, Pc 

97 


1692-1706 

15 

I 

Winteraceae/Canellaceae 

98 


1692-1695 

4 

u 

Tl 

99 


1700-1701 

2 

u 

Tp 

l(K) 


1707-1722 

16 

I 

Winteraceae/Canel laceae 

101 


1709 

1 

u 

Ce 

102 


1713-1716 

4 

u 

Cw 

1 ()3 


1717-1718 

2 

H 

Cw, Pm 

104 


1719-1720 

2 

U 

Ce 

105 


1726-1727 

2 

U 

Tl 

106 


1730 

1 

1 

W interaceae/Canellaceae 

107 


1737 

1 

U 

Pc 

108 


1738 

1 

U 

Pa 

109 


1736-1748 

10 

u 

Pc 

110 


1750-1759 

11 

1 

Winteraceae/Canellaceae 

111 


1757-1758 

2 

u 

Tp 

112 


1763-1764 

2 

H 

Cw, Cd, Pm 

113 


1765-1774 

10 

I 

Cd, Cm, Ws 

114 


1775 

1 

H 

Cd, cm 

115 


1776-1778 

3 

U 

Cd 


2 (bp 1690—1722), the alignment between families 
was also relatively straightforward. The two prob¬ 
lematic regions were partitioned equally between 
the two families as non-overlapping separate indels 
(e.g., positions 1143—1153 were coded as gaps in 
Winteraceae and positions 1154—1163 were coded 
as gaps in Canellaceae). The effect of this was to 
reduce the probability of incorrect homology as¬ 


sessment between families while retaining potential 
phylogenetic information within each family. 

Putative secondary structures for ITS 2 were de¬ 
termined for each sequence, and these structures 
were generally consistent with the final alignment 
(data not shown). An apparent contradiction be¬ 
tween secondary structure and DNA alignment was 
found in the v5 region (see Hershkovitz & Zimmer 
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[1996] for details) and corresponds to bp 1704— 
1734. Although this region consistently formed a 
bulge in secondary structure analyses across all in¬ 
dividuals, the base composition and sequence 
length clearly differed between families (11—15 bp 
with 46—60% G/C in Winteraceae and 7—16 bp 
with 75—85% G/C in Canellaceae). This result fur¬ 
ther supports treatment of this region as non-over¬ 
lapping indel events. No known secondary structure 
exists for ITS 1; consequently, structural evidence 
supporting the alignment of ITS 1 could not be gen¬ 
erated. 

A total of 115 indels were identified in the ITS 
alignment, of which 53 (46%) were parsimony-in¬ 
formative, 38 (33%) were autapomorphic, and the 
remaining 24 (21%) were homoplastic. Table 2 pre¬ 
sents these indels along with position information 
(character number in the data set), indel length, 
and the species where indels were found. 

Chloroplast ImL-F. The frnL-F spacer regions 
were easily aligned by eye within and between fam¬ 
ilies. Overall, a total of 24 indels (Table 2: A—X) 
were identified in the alignment. Of them, 9 
(37.5%) were parsimony-informative, 12 (50.0%) 
were autapomorphic, and the remaining 3 (12.5%) 
were homoplastic. Especially noteworthy was indel 
P (Table 2) with 105 bases present in the Winter¬ 
aceae, including Takhtajania, but this indel was not 
present in the Canellaceae. 

Combined data set. Two partitions of the com¬ 
bined alignment were defined, (1) trnC-Y (bp 1— 

987) and (2) ITS (bp 988-1778). These partitions 

were tested for congruence with the data partition- 
homogeneity test (1000 replicates) and yielded a P 
value of 0.367. This result fails to reject the null 
hypothesis of congruence between the two parti¬ 
tions. Therefore, all phylogenetic analyses were 
performed on both the individual and combined 
data sets. 

PATTERNS OF SEQUENCE EVOLUTION 

The aligned combined data set was 1778 char¬ 
acters long. Positions 1 to 987 represented the 
trnL-F chloroplast sequence region; positions 988 
to 1778 represented the ITS nuclear ribosomal se¬ 
quence region including the 5.8S coding region. 
Unlike the situation for certain genera of Wintera¬ 
ceae (Belliolum, Bubbia, Exospermum, and Zygo¬ 
gynum), no polymorphisms were detected in the 
ITS sequence for Takhtajania or any member of the 
Canellaceae. 

The ITS 1 sequences in Winteraceae ranged 
from 238 bp in Takhtajania to 252 bp in Drimys. 
In Canellaceae, ITS 1 sequences ranged from 244 


bp in Capsicodendron to 274 bp in Pleodendron. 
The 5.8S coding region was 164 bp in all individ¬ 
uals of both families, although it was only 163 bp 
long in Canella. ITS 2 sequences in Winteraceae 
ranged from 215 bp in Pseudowintera colorata to 
228 bp in Belliolum, with Takhtajania having a 222 
bp sequence. Canellaceae ITS 2 sequences ranged 
from 187 bp in Cinnamodendron to 211 bp in Ca¬ 
nella. 

The Winteraceae chloroplast group I intron se¬ 
quences ranged from 490 bp in Tasmannia insipida 
to 508 bp in Takhtajania perrieri. In Canellaceae, 
these sequences ranged from 477 bp in Canella 
winterana to 495 bp in Cinnamosma madagascar- 
iensis. The 5' trnh exon sequences were 50 bp in 
all individuals of both families. The intergenic 
spacer sequences in Winteraceae ranged from 360 
bp in Drimys to 368 bp in both Tasmannia insipida 
and Takhtajania perrieri. All Canellaceae sequenc¬ 
es were markedly shorter, ranging from 247 bp in 
Warburgia salutaris to 259 bp in Canella winter- 
ana, Cinnamodendron ekmanii, and Pleodendron 
macranthum. The 40 bp partial trnF exon sequenc¬ 
es were identical in length in all individuals of both 
families. 

PHYLOGENETIC ANALYSES 

Maximum parsimony 

Nuclear ITS. Using the Branch and Bound 
search procedure in PAUP*, with equal weighting 
of positions and Canellaceae selected as the out¬ 
group, a single most parsimonious tree of 489 steps 
(366 steps excluding uninformative characters) was 
obtained with the ITS data set. The consistency in¬ 
dex (Kluge & Farris, 1969) was 0.8221 (0.7623 

excluding uninformative characters) and retention 
index (Farris, 1989) was 0.8793. Overall, 317 var¬ 
iable nucleotide positions were counted; of these, 

206 (or 26% of the 791 total characters) provided 
parsimony-informative changes. 

To assess robustness of these results, the bootstrap 
procedure (1000 replicates) and decay analyses were 

run in PAUP* (TOPOLOGICAL CONSTRAINTS 

were used to determine decay values for each node). 
These results are shown in Figure 1. 

A polytomy including Belliolum, Exospermum, 
and Zygogynum (Zygogynum s.l.) was supported as 
monophyletic and sister to Bubbia by a 97% boot¬ 
strap value and decay index of + 4. The Bubbia + 
Zygogynum s.l. clade was monophyletic and sister 
to Pseudowintera (represented by two species, 
monophyletic in 99% of the bootstrap iterations 
[decay = + 5]), supported by a 100% bootstrap 
value (decay = + 10). Drimys was sister to this 
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95/99 
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1/1 

1 63/62 
♦ 1 M 

1/1 





64/62 

♦ 1/+1 


5/6 

4/4 




97/99 


1/1 


+ 4/+4 

3/3 


5/7 


2/3 


f-H 



24/31 


9/10 



28/40 


21/23 




7/8 


21/24 


26/27 



84/85 

♦ 4/+4 


88/93 

*71*7 


73/77 

+ 3/+3 


33/39 


I I I I l 

II 111 


34/39 


5/5 



35/37 


100/100 
* 251*27 


9/10 



51/56 


Zygogynum bicolor 


Belliolum pancheri 


Zygogynum balansae 


Exospermum stipitatum 


Zygogynum acsmithii 


Bubbia comptonii 


5/5 


99/100 

♦ 5/+5 1 

3/3 


Pseudowintera axillaris 


Pseudowintera colorata 


Drimys winteri 


95/128 


, 11/18 


100/100 

1 99/100 

7/9 

+ 61/+91 

+ 8/+13 



Tasmannia insipida 


Tasmannia lanceolata 


Takhtajania perrieri 



Warburgia salutaris 


Cinnamosma madagascariensis 


Capsicodendron denisii 



Cinnamodendron ekmanii 



Pleodendron macranthum 


Can el la winterana 


figure 1. Single most parsimonious tree generated both from the I TS data set alone and the combined ITS//rnl,-F 
data set. Tree length for IT'S tree = 489 steps (366 steps excluding uninformative characters), consistency index = 
0.8221 (0.7623 excluding uninformative characters), and retention index = 0.8793. Length of combined tree = 593 
steps (434 steps excluding uninformative characters), consistency index = 0.8381 (0.7788 excluding uninformative 
characters), and retention index = 0.8994. Numbers to the left of the slash represent those generated from ITS data 
alone, while numbers to the right of the slash represent those generated from the combined data set. Branch lengths 
shown above branches, bootstrap values (1000 replicates) and decay values shown below branches, respectively. Dark 
bars represent non-homoplastic insertion/deletion events inferred from the ITS alignment (see Table 2). 
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assemblage in 95% of the bootstrap iterations (de¬ 
cay = + 3). Tasmannia (represented by two spe¬ 
cies) was monophyletic with 99% bootstrap support 
(decay = + 8) and sister to the Drimys + Pseu- 
dowintera + Bubbia + Zygogynum s.l. clade in 
86% of the bootstrap iterations (decay = + 3). 
Takhtajania was found sister to the remainder of 
the Winteraceae with a bootstrap value of 94% (de¬ 
cay = + 8). 

In the Canellaceae, Cinnamosma and Warburgia 
formed a clade sister to Capsicodendron in 73% of 
the bootstrap replicates, with a decay index of + 
3. Cinnamodendron and Pleodendron formed a 
well-supported clade (100% of the bootstrap rep¬ 
licates, decay = + 25) sister to the Capsicodendron 
+ Cinnamosma + Warburgia clade in 88% of the 
replicates (decay = + 7). Finally, Canella was 
found sister to the remainder of the family with 
84% bootstrap support (decay = + 4). 

Chloroplast trnL-F. Using the Branch and 
Bound search procedure in PAUP*, with equal 
weighting of positions and Canellaceae selected as 
the outgroup, two most parsimonious trees of 103 
steps (67 steps excluding uninformative characters) 
were obtained with the trnC-Y data set. The con¬ 
sistency index was 0.9223 (0.8806 excluding un¬ 
informative characters), and the retention index was 
0.9657. Overall, 90 variable nucleotide positions 
were counted; of these, 55 (or 5.6% of the 987 total 

characters) provided parsimony-informative chang¬ 
es. These two trees are presented in Figure 2 along 
with bootstrap values and decay indices. 

With the exception of alternative equally parsi¬ 
monious basal relationships, the topology generated 
by the trnY-Y sequence data for Winteraceae was 
not in conflict with the tree generated by the ITS 
data, albeit there were both lower resolution and 
lower bootstrap support values. The polytomy of 
Belliolum, Exospermum, and Zygogynum ( Zygo¬ 
gynum s.l.) expanded to include Bubbia (78% boot¬ 
strap and decay = +2). The two species of Pseu- 
dowintera were unresolved with respect to each 
other and to the Bubbia + Zygogynum s.l. poly¬ 
tomy, and together these were only weakly sup¬ 
ported as monophyletic and sister to Drimys (69% 
bootstrap and decay = + 1). The two species of 
Tasmannia were strongly supported as sister taxa 
in 98% of the bootstrap iterations (decay = + 5). 
Finally, the frnL-F sequence data could not resolve 
basal relationships within the family; both Takhta¬ 
jania (Fig. 2a) and Tasmannia (Fig. 2b) as sister to 
the remaining Winteraceae were equally parsimo¬ 
nious. 

Within Canellaceae, except for a strong sister re¬ 
lationship between Cinnamodendron and Pleoden¬ 


dron (88% bootstrap and decay = + 2), the topol¬ 
ogy generated by frnL-F data resulted in much 
lower levels of bootstrap support, and, moreover, 
suggested a fundamental split between Old World 
and New World lineages. Thus, rather than being 
sister to the remaining Canellaceae as suggested by 
the ITS-based phylogeny, Canella was nested with¬ 
in a weakly supported clade of New World genera. 
To force such a division of the family into Old 
World and New World clades requires an additional 
11 steps in the ITS-based phylogeny. 

Combined data set. Using the Branch and 
Bound search procedure in PAUP*, with equal 
weighting of positions and all Canellaceae genera 
selected as outgroups, a single most parsimonious 
tree of 593 steps (434 steps excluding uninforma¬ 
tive characters) was obtained with the ITSArnL-F 
combined data set. The consistency index was 

0.8381 (0.7788 excluding uninformative charac¬ 
ters) and the retention index was 0.8994. The to¬ 
pology derived from this combined data set was 
identical to that generated from the ITS data alone. 
The branch lengths, bootstrap values, and decay 
indices are presented in Figure 1. 

Overall, with the addition of the trnL-F data, 
bootstrap values either remained similar or in¬ 
creased in comparison with the ITS-based phylog¬ 
eny. Within the Winteraceae, the most notable in¬ 
creases resulted in stronger support for basal 
relationships in the family. The position of Takh¬ 
tajania as basal and sister to the rest of the family 
was further strengthened (ITS = 94% bootstrap and 
+ 8 decay index vs. combined data = 99% boot¬ 
strap and +11 decay index), as was the position 
of Tasmannia as the next branch within the family. 
The combined data set provided especially strong 
support for the inclusion of Drimys in the clade 
consisting of all other members of the family and 
sister to Tasmannia (ITS = 86% bootstrap and + 
3 decay index vs. combined data = 97% bootstrap 
and + 6 decay index). Within the Canellaceae, 
bootstrap support also remained similar or in¬ 
creased only slightly with the addition of the trnL- 

F data (Fig. 1). 

Maximum likelihood 

Nuclear ITS. For the ITS data set, the best 
model of DNA substitution was the general time- 
reversible model with a proportion of sites assumed 
to be invariable and equal rates assumed across 
variable sites (GTR + I). With this model, param¬ 
eters were fully optimized using the tree generated 
from the parsimony search (Fig. 1). With the model 
parameters defined, a heuristic search was con- 
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Figure 2a, b. Two most parsimonious trees generated using trnL-V data. Tree length = 103 steps (67 steps excluding 
uninformative characters), consistency index = 0.9223 (0.8806 excluding uninformative characters), and retention index 
= 0.9657. Branch lengths shown above branches, bootstrap values (1000 replicates) and decay values shown below 
branches, respectively. Dark bars (Fig. 2a) represent non-homoplastic insertion/deletion events inferred from the trnh- 
V alignment (see Table 2). 
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Figure 2. Continued. 
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Figure 3a. One of four most likely phylograrns generated with ITS data ( — In = 3371.16) with bootstrap values 
(1(HH) replicates) above branches. 
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Figure 3b. Single most likely phylogram generated with trnL-V data ( — In = 1996.3473) with bootstrap values 
(1(X)0 replicates) above branches. Figure 3c is on page 428. 


ducted that resulted in three trees, one of which 
was identical in topology to the most parsimonious 
tree (Fig. 1). These three trees were used to further 
optimize model parameters. The resulting parame¬ 
ter values did not differ from those generated by 
using the single maximum parsimony tree. Thus, 


ten heuristic searches (with random taxon addition, 
multrees, and steepest descent active) were con¬ 
ducted using the parameters defined in the initial 
search. These ten random taxon addition searc hes 
generated a total of four equally likely trees (—In = 
3371.16), one of which is presented in Figure 3a. 



























428 


Annals of the 

Missouri Botanical Garden 



Figure 3c. Single most likely phylogram generated with combined ITS//raL-F data set (In 
bootstrap values (1000 replicates) above branches. 


5508.26165) with 


In all four trees the relationships within Winter- 
aceae were identical to those generated by the MP 
search, with Takhtajania sister to the remainder of 
the Winteraceae. The four trees differed in the po¬ 
sition of Canella in the Canellaceae. Either Canella 
was sister to the Warburgia + Cinnamosma + Cap - 
sicodendron clade, or sister to the Cinnamodendron 
+ Pleodendron clade, or sister to the remainder of 


the family, or, finally, in an unresolved polytomy 
with the three above-mentioned clades. Although 
Canella is included within a monophyletic Canel¬ 
laceae in 100% of the bootstrap iterations, boot¬ 
strap values do not support a resolved phylogenetic 
position of Canella within the family. 

Chloroplast tm L-F. For the trnh-F data set, the 
best model of DNA substitution was the HKY85 
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model with a proportion of sites assumed to be in¬ 
variable and equal rates assumed across variable 
sites (HKY85 + I). With this model, parameters 
were fully optimized using the two trees generated 
from the parsimony search (Fig. 2a, b). With the 
model parameters defined, a heuristic search was 
conducted that resulted in a single tree identical in 
topology to one of the four most likely trees gen¬ 
erated from the maximum likelihood ITS search 
(Fig. 3a). Using this tree, parameters were re-esti- 
mated and the resulting optimized values were used 
for ten heuristic searches (random taxon addition, 
multrees, and steepest descent active). The result 
of these searches generated a single most likely tree 

( — In = 1996.3473), which is presented in Figure 

3b along with bootstrap values for each node. 

The single tree from the maximum likelihood 
search was identical to one of the two trees gen¬ 
erated by the maximum parsimony analysis using 
trnL-F with one exception. The placement of Ca- 
nella sister to Pleodendron and Cinnamodendron 
was not supported. Rather, Capsicodendron , Cin- 
namosma , and Warburgia were supported as an un¬ 
resolved trichotomy in 69% of the bootstrap itera¬ 
tions. Cinnamodendron and Pleodendron were 
supported as sister taxa in 87% of the bootstrap 
iterations, and finally Canella was found in an un¬ 
resolved trichotomy with the two above-mentioned 
clades in a monophyletic Canellaceae in 100% of 
the bootstrap replications. Within the Winteraceae 
an unresolved polytomy including Belliolum , Bub - 
bia , Exospermum , and Zygogynum was found to be 
sister to the two species of Pseudowinter a in 71% 
of the bootstrap replicates. In 78% of the bootstrap 
iterations, Drimys joined the above in a monophy¬ 
letic clade, but was only weakly supported by a 
54% bootstrap value as sister to them. There was 
less than 50% bootstrap support for the basal po¬ 
sition of Tahhtajania sister to the remainder of Win¬ 
teraceae. 

Combined data set . For the combined data set, 
the best model of DNA substitution was the general 
time-reversible model. Among site rate variation 
was accommodated for by dividing the data into two 
partitions (ITS and trnL-F), thus allowing among 
site rate variation to be estimated for each molecule 
separately (GTR + SS). With this model, parame¬ 
ters were optimized using the tree generated from 
the maximum parsimony search (Fig. 1). With the 
model parameters defined, a heuristic search was 
conducted that resulted in a single most likely tree 
identical in topology to the most parsimonious tree 
using the same data set (Fig. 1). This tree was used 
to further optimize model parameters. The resulting 
parameter values did not differ, and thus ten heu¬ 


ristic searches (random taxon addition, multrees, 
and steepest descent active) were conducted using 
the parameters defined in the initial search. These 
ten random taxon addition searches generated a 

single most likely tree (In = —5508.26165). The 
bootstrap values are presented in Figure 3c. 

In general, with the exception of support for Ca¬ 
nella as basal within Canellaceae, the bootstrap 
values were consistent with those generated by the 
parsimony analyses of the combined data set. In 
the parsimony analysis, Canella was sister to the 
remainder of the Canellaceae in 85% of bootstrap 
iterations, whereas the likelihood analyses recov¬ 
ered this relationship in only 53% of the bootstrap 
replicates. 

Discussion 

In the initial molecular phylogenetic study of 
Winteraceae, Suh et al. (1993) found it difficult to 
align the ITS region with certainty for several ex¬ 
emplar outgroup species. The present study in¬ 
creased generic sampling in Canellaceae from a 
single exemplar ( Canella ) to a representative of all 
genera, which greatly facilitated homology assess¬ 
ment between the two families. In addition, the 
more slowly evolving trnL-F chloroplast regions, for 
which no difficulty in alignment was encountered, 
were investigated. With Canellaceae as the out¬ 
group, the rooting of Winteraceae is no longer am¬ 
biguous. Phylogenetic analyses of ITS and com¬ 
bined ITS/trnL-F DNA sequences provided a 
well-resolved molecular phylogeny for the Winter¬ 
aceae that placed Tahhtajania sister to the remain¬ 
ing genera in the family. Analysis of trnL-F se¬ 
quence data alone failed to resolve the basal 
topology within the family, with either Tahhtajania 
or Tasmannia equally parsimonious as the basal 
branch. Nevertheless, the addition of trnL-F se¬ 
quences to a combined ITS/trnL-F analysis resulted 
in increased support values for the same topology 
derived from ITS sequences alone. Maximum like¬ 
lihood analyses of each of the three data sets mir¬ 
rored those obtained through maximum parsimony. 

ITS and trnL-F sequences provided generally 
weaker support for phylogenetic relationships with¬ 
in Canellaceae. Although parsimony analysis of ITS 
data revealed moderate support for a basal position 
of Canella sister to the remaining members of the 
family, parsimony analysis of trnL-F data weakly 
suggested a fundamental New World/Old World 
split in Canellaceae. Parsimony analysis of the 
combined molecular data set resulted in a topology 
identical to that generated by ITS data alone, with 
only slight increases in bootstrap support and iden- 
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tical decay values. Maximum likelihood analyses of 
each of the three data sets further demonstrated the 
poor resolution these molecules provide for esti¬ 
mating basal relationships in Canellaceae. Never¬ 
theless, ITS, combined ITS/JrnL-F, and all maxi¬ 
mum likelihood analyses strongly suggest that the 
southeastern Brazilian endemic Capsicodendrort 
forms a clade with the African/Malagasy genera 
Warburgia and Cinnamosma. 

The alignment of ITS and fmL-F sequences for 
the 12 species of Winteraceae and 6 species of Ca¬ 
ne llac •eae revealed numerous probable insertion- 
deletion events. When mapped on the trees, these 
indels furnished additional support for certain 

clades. A total of 37 indels (29 in ITS [Fig. 1] and 
8 in trnL-F [Fig. 2a]) support the separation of Win¬ 
teraceae from Canellaceae. Three ITS indels sup¬ 
port the monophyly of Pseudowinter a, and three dif¬ 
ferent ITS indels support the monophyly of the 
Zygogynum s.l. (Bubbin, Belliolum , Exospermum , 
and Zygogynum) assemblage. One indel each in 
both nuclear and chloroplast sequences (Figs. 1 
and 2a, respectively) support inclusion of Drimys 
in the Zygogynum s.l. + Pseudowinter a clade, and 
thus provide additional evidence of its distinctness 
from Tasmannia. Within Canellaceae, four ITS in¬ 
dels support the sister relationship between Cin- 
namodendron and Pleodendron , whereas a single 
ITS indel supports a sister relationship between Af¬ 
rican Warburgia and Malagasy Cinnamosma. The 
basal position of Canella in the family is further 
suggested by three ITS indels that support the 
monophyly of the remaining genera of Canellaceae. 
None of the indels found in the trn L-F Canellaceae 
sequences were phylogenetically informative. In to¬ 
tal, 139 indels were scored overall and 114 (82%) 
were consistent with the trees generated from the 
sequence data. The remaining 25 were homoplastic 
and in some cases difficult to score. For example, 
in trnL-F sequences, long adenine runs occurred at 
varying numbers in closely related taxa. The length 
of these runs was not consistent when mapped on 
the tree, and delimiting homology was difficult 
when scoring this region. No indels were detected 
in either ITS or trnL-F that bear on the position of 
Takhtajania as sister to the remainder of Wintera¬ 
ceae. 

Zygogynum s.l.. Pseudowintera, and Takhtajania 
exhibit a short, early-rupturing involucrum (the 
congenitally fused two, or rarely three, outermost 
perianth parts, which could as well be considered 
the “calyx”) that is persistent in fruit vs. the con¬ 
dition in both Drimys and Tasmannia where the 
ultimately caducous, late-rupturing involucrum be¬ 
comes as long as the fully developed next inner 


tepal whorl, and therefore completely encloses the 
flower just prior to anthesis (Doust, 2000 this is¬ 
sue). A persistent, short calyx occurs in Canella¬ 
ceae that is similar to the ruptured, persistent in¬ 
volucrum of basal Takhtajania , suggesting that 
Pseudowintera and Zygogynum s.l. have retained 
the ancestral plesiomorphic condition of an early- 
rupturing involucrum. In conjunction with the hy¬ 
pothesized phylogenetic positions of Drimys and 
Tasmannia , the differing ontogenetic patterns of te¬ 
pal initiation and development exhibited by these 
two genera (Doust, 2000) imply that the late-rup¬ 
turing, caducous involucrum has evolved twice in¬ 
dependently. Similarly, deep red tepal color, pre¬ 
sent in Canella and other Canellaceae, as well as 
Takhtajania and many Zygogynum s.l. species, 
may be the ancestral condition, with the evolution 
of white tepals through the loss of pigmentation 
and/or the acquisition of an optical tapetum (En- 
dress et al., 2000 this issue) occurring a number of 
times independently. 

The historical biogeographical implications of 
Takhtajania in a basal position sister to the re¬ 
maining genera of the Winteraceae have been dis¬ 
cussed in the context of the distribution of fossil 
Winteraceae (Doyle, 2000 this issue) and ecophys- 
iological constraints (Feild et al., 2000 this issue). 
It is reasonable to assume that the basal branch 
leading to Takhtajania became isolated in Mada¬ 
gascar after reaching there from continental Africa, 
and also that Madagascar may have played a prom¬ 
inent role in the migration of the family to Aus¬ 
tralasia and South America via Antarctica during 
the mid-Cretaceous. The presence of two different 
winteraceous pollen types in the Miocene of the 
southwestern Cape possibly referable to Tasmannia 
and the most advanced Zygogynum s.l. clade is 
more puzzling and problematic (Coetzee & Prag- 
lowski, 1988). With South Africa relatively isolated 
from both South America and Antarctica from the 
mid-Cretaceous onward (Smith et al., 1994), either 
all of the currently extant clades of Winteraceae 
had already evolved prior to that point in time, or 
Madagascar may have continued to serve as a con¬ 
duit for migration of more advanced Winteraceae 
back into South Africa from Australasia. In com¬ 
parison to the Southern Gondwanan pattern of ra¬ 
diation in Winteraceae, the preliminary phyloge¬ 
netic hypothesis of infrafamilial relationships 
within Canellaceae suggests an “inverted” Northern 
Gondwanan biogeographical history. Basal clades of 
Canellaceae are centered in northern South Amer¬ 
ica and the West Indies, with a more advanced 
clade exhibiting a southeastern Brazil/African—Mal¬ 
agasy split related to the middle Cretaceous sepa- 
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ration of South America and Africa (Goldblatt, 
1993) and a subsequent dispersal event to Mada¬ 
gascar. Thus, the modern-day sympatric occurrence 
of the two Malagasy endemic genera Cinnamosma 
(Canellaceae) and Takhtajania (Winteraceae) may 
well be the end result of vastly divergent biogeo¬ 
graphic histories. 
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